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Abstract: We study the 6 dependence of the spectrum of four-dimensional SU(iV) 
gauge theories, where 9 is the coefficient of the topological term in the Lagrangian, 
for N > 3 and in the large- iV limit. We compute the 0{9 2 ) terms of the expansions 
around 9 = of the string tension and the lowest glueball mass, respectively a{9) = 
a (1 + s 2 9 2 + ...) and M(9) = M (1 + g 2 9 2 + ...), where a and M are the values 
at 9 = 0. For this purpose we use numerical simulations of the Wilson lattice 
formulation of SU(iV) gauge theories for N = 3,4, 6. The 0{9 2 ) coefficients turn out 
to be very small for all iV > 3. For example, s 2 = —0.08(1) and g 2 = —0.06(2) for 
N = 3. Their absolute values decrease with increasing N. Our results are suggestive 
of a scenario in which the 9 dependence in the string and glueball spectrum vanishes 
in the large- N limit, at least for sufficiently small values of \9\. They support the 
general large- N scaling arguments that indicate 9 = 9/N as the relevant Lagrangian 
parameter in the large- N expansion. 
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1. Introduction 

Four- dimensional SU(iV) gauge theories have a nontrivial dependence on the angle 
9 that appears in the Euclidean Lagrangian as 

C e = \f; v {x)F; v {x) - i9^e, vpa F^(x)F^(x) (1.1) 

(q(x) = ■^s^fivpaF^ l/ (x)Fp (T (x) is the topological charge density). Indeed, the most 
plausible explanation of how the solution of the so-called U(1)a problem can be 
compatible with the 1/N expansion (performed keeping g 2 N fixed |IJ) requires a 
nontrivial 9 dependence of the ground-state energy density F(9), 

F{9) = -i In J [dA] exp {- j d d xC^ , (1.2) 

in the <i-dimensional pure gauge theory to leading order in 1/N @, [J. The large- N 
ground-state energy is expected to behave as [||, |. [| 

AF{6) = F{9) - F(0) =A9 2 + (l/N 2 ) (1.3) 

for sufficiently small values of 9, i.e. 9 < ir. This has been supported by Monte 
Carlo simulations of the lattice formulation of SU(iV) gauge theories 0. Indeed, the 
numerical results for N = 3, 4, 6 are consistent with a scaling behavior around 9 = 
given by 

f{9) = a- 2 AF(9) = l -C9 2 (l + b 2 9 2 + ...), (1.4) 
C = Co, + c 2 /N 2 + b 2 = b 2 , 2 /N 2 + ... 
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where a is the string tension at 9 — 0. C is the ratio x/ a2 where \ = J d 4 x(q(x)q(0)) 
is the topological susceptibility at 9 = 0. Its large- N limit Coo is 0, « 0.022. 
Moreover, estimates of c 2 and 6 2)2 are c 2 ~ 0.06, and 6 2j2 ~ —0.2 {b 2 ~ —0.02 for 
SU(3) [0, |). Note that Eq. (|ll]) can be recast in the form 



f(9) = N 2 f(9 = 9/N), (1.5) 
]{9) = \C9\1 + W + ...), 

where 6 2 = 6 2j2 + 0(1/N 2 ) = 0(1). This is consistent with general large- N scaling 
arguments applied to the Lagrangian (|1 . 1| ) , which indicate 9 = 9/N as the relevant 
Lagrangian parameter in the large- N limit of the ground-state energy |3J. 

Another interesting issue concerns the 9 dependence of the spectrum of the the- 
ory. This is particularly interesting in the large- N limit where the issue may also be 
addressed by other approaches, such as AdS/CFT correspondence applied to non- 
supersymmetric and non conformal theories, see e.g. Ref. |TI|. The analysis of the 



9 dependence of the glueball spectrum using AdS/CFT suggests that the only effect 
of the 9 term in the leading large- A limit is that the lowest spin-zero glueball state 
becomes a mixed state of ++ and h glueballs, as a consequence of the fact that 



the 9 term breaks parity, but its mass does not change ||11|| . 

In this paper we present an exploratory numerical study of the 9 dependence in 
the spectrum of SU(A) gauge theories. For this purpose we use numerical simulations 
of the Wilson lattice formulation. Numerical Monte Carlo studies of the 9 dependence 
are made very difficult by the complex nature of the 9 term. In fact the lattice action 
corresponding to the Lagrangian (|1.1| ) cannot be directly simulated for 9^0. Here 
we restrict ourselves to the region of relatively small 9 values, where one may expand 
the observable values around 9 = 0. We consider the string tension and the lowest 
glueball mass. We write 

a(9) = a (l + s 2 9 2 + ...), (1.6) 

where a is the string tension at 9 = 0. When N > 4 analogous expressions can be 
written for the other independent fc-strings associated with group representations of 
higher n-ality. Moreover, for the lowest glueball state we write 

M(9) = M(l + g 2 9 2 + ...) (1.7) 

where M is the ++ glueball mass at 9 = 0. Then the coefficients of these expansions 
can be computed from appropriate correlators at 9 = 0. The 0(9 2 ) coefficients s 2 and 
<? 2 are dimensionless quantities, which should approach a constant in the continuum 
limit, with 0(a 2 ) scaling corrections. The idea is analogous to the one exploited in 
Ref. to study the ^-dependence of the ground-state energy. 

We shall present results for four- dimensional SU(A) gauge theories with N = 
3,4,6. The estimates of the 0(9 2 ) coefficients turn out to be very small for all 
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N > 3. For example S2 = —0.08(1) and gi = —0.06(2) for N = 3. Moreover, their 
absolute values decrease with increasing N. We also observe that the 0(9 2 ) terms 
are substantially smaller in dimensionless ratios such as M/y/a and, for iV > 3, the 
ratios of independent k strings, = cy k /a. Our results are suggestive of a large- N 
scenario in which the 9 dependence in the string and glueball spectrum vanishes 
around 9 = 0. They are consistent with the general large- N scaling arguments 
indicating 9 = 9/N as the relevant parameter in the large- N limit. We also show 
that a similar scenario emerges in the two-dimensional CP^ -1 models by an analysis 
of their 1/N expansion. 

The paper is organized as follows. In Sec. @ we outline the numerical method 
to estimate the 0(6 2 ) terms of the expansion in powers of 9 around 9 = 0. The 
results of our exploratory numerical study are presented in Sec. [| Finally, in Sec. [| 
we discuss the 9 dependence of two-dimensional CP^ -1 models within their 1/N 
expansion around their large- N saddle-point solution. 



2. Numerical method 
2.1 Monte Carlo simulations 

We consider the Wilson formulation of lattice gauge theories: 

S = -N/3 Tr [U^x)U u {x + /i)Ul(x + v)Ut(x) + h.c] , (2.1) 

X,fl>V 

where U^(x) G SU(iV) are link variables. In our simulations we employed the 
Cabibbo-Marinari algorithm |l2j to upgrade SU(iV) matrices by updating their SU(2) 
subgroups (we selected iV(iV — l)/2 subgroups and each matrix upgrading consists 
of iV(iV — l)/2 SU(2) updatings). This was done by alternating microcanonical 
over- relaxation and heat-bath steps, typically in a 4:1 ratio. 

Computing quantities related to topology using lattice simulation techniques is 
not a simple task. In a lattice theory the fields are defined on a discretized set, 
therefore the topological properties are strictly trivial. One relies on the fact that 
the physical topological properties are recovered in the continuum limit. Various 
techniques have been proposed and employed to associate a topological charge Q 
to a lattice configuration, see, e.g, Refs. O, 114] for techniques based on bosonic 



operators, and Refs. [15|, [16[ for techniques based on fermionic estimators. The 
most robust definition of topology on the lattice is the one obtained using the index 
of the overlap Dirac operator. However, due to the computational cost of fermionic 
methods and the need for very large statistics to measure correlations of Polyakov and 
plaquette operators with topological quantities, we decided to use the simpler cooling 
method, implemented as in Ref. 0. Direct comparison with a fermionic estimator 
is known to show a good agreement in the case of SU(3) ^Hh 0, |T7[, supporting the 
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idea that the cooling method is fairly stable in this case. Moreover, the agreement 
among different methods is expected to improve with increasing iV fll8l , |i~9| . 

A severe form of critical slowing down affects the measurement of Q, posing a 
serious limitation for numerical studies of the topological properties in the continuum 
limit, especially at large values of N. The autocorrelation time tq of the topological 
modes rapidly increases with the length scale, much faster than the standard square 
law of random walks ]7|, ^D|. The available estimates of tq appear to increase as 
an exponential of the length scale, or with large power laws. A qualitative expla- 
nation of this severe form of critical slowing down may be that topological modes 
give rise to sizeable free-energy barriers separating different regions of the configu- 
ration space. As a consequence, the evolution in configuration space may present a 
long-time relaxation due to transitions between different topological charge sectors. 
This dramatic effect has not been observed in plaquette-plaquette or Polyakov line 
correlations, suggesting an approximate decoupling between topological modes and 
nontopological ones, such as those determining the confining properties and the glue- 
ball spectrum. But, as we shall see, such a decoupling is not complete. Therefore 
the strong critical slowing down that is clearly observed in the topological sector will 
eventually affect also the measurements of nontopological quantities, such as those 
related to the string and glueball spectrum. 

2.2 The 0(9 2 ) coefficients of the 9 expansion 

In this subsection we describe the method to determine the 0(8 2 ) coefficients of 
the 9 expansions such as Eqs. ( |1.6|) and (|1.7|) . Let us first discuss the case of the 
fundamental string tension. The string tension can be determined from the torelon 
mass, i.e. the mass describing the large-time exponential decay of the wall-wall 



correlations Gp of Polyakov lines [21|. In the presence of a 9 term 



G P (t, 9) - (A P {t))g - f [du]e „j d i xCe (2-2) 

where 

A P (t) = TrP(0;0) Tr P{x x , x 2 ; t), (2.3) 

X\,X2 

and P(xi,X2',t) is the Polyakov line along the x^ direction of size L. The time sep- 
aration t is an integer multiple of the lattice spacing a: t = n t a. The correlation 
Gp can be expanded in powers of 9. Here, we are considering the case of the funda- 
mental string tension, but the discussion can be easily extended to any other group 
representation, by replacing the trace with the corresponding character in Eq. ( |2.3| ). 
Taking into account the parity symmetry at 9 = 0, we obtain 

G P (t,9) = G { °\t) + h) 2 Gf{t) +0(9*), (2.4) 
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where 



G$\t) = (A P (t)) e=0 , 
cf(t) = -(A P (t)Q 2 )e=o+ (A P (t)) e=0 (Q 2 ) e=0 



[2.5) 
(2.6) 



and Q is the topological charge. 

The correlation function Gp is expected to have a large-t exponential behavior 



G P (t,9) ^ B(e)e- E{e)t , 



(2.7) 



where E(6) is the ^-dependent energy of the lowest state (torelon mass), and B(6) 
is the overlap of the source with the lowest-energy state. If the lattice size L is 
sufficiently large, the lowest-energy states describing the long-distance behavior of 
Polyakov correlators should be those of a string-like spectrum. Then, the string 
tension is extracted using the relation 



E(9) = a(9)L 



71 

3L 



(2.8) 



Here we are assuming that the 0(1/ L) (Liischer) correction is independent of 9. 
Actually we are also assuming the so called free string spectrum 



a 



LI 



7i 

3/1 



71- 



l a = y/aL 



(2.9) 



(n: excitation level), obtained neglecting the self- interaction terms in the string ef- 
fective action, see e.g. Ref. |]2"2"|| . As shown in Ref. j2"2"fl , only the 0(1/ L) correction 
should be universal, while subleading corrections are generally expected. They de- 
pend on the unknown coefficients of the higher order terms of the effective QCD 
string action. For example, besides the free string spectrum, one may also consider 



the Nambu-Goto string spectrum |22| , |23 

W n = oL |1 



2tt 8tt 

m + "I 



1/2 



(2.10) 



In particular, if one assumes the Nambu-Goto spectrum, instead of Eq. ( |2.8| ), one 
should use the Nambu-Goto lowest-energy state to determine the string tension, i.e. 



Wn 



2tt 

" L 1 1 " 51 



1/2 



oh 



(2.11) 



Therefore the Nambu-Goto string spectrum leads to a different estimate of the string 
tension at finite l a : 



Cfs — o"ng 



1 - 



7T" 



18/! 



0(V 



(2.12) 
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where Of s and cxng are the string tensions extracted from the torelon mass assuming 
respectively the free and Nambu-Goto spectrum. 

Lattice sizes such that l a > 3 should be sufficiently large to have an effective 
string picture parametrized by a constant string tension a [p4" |. Then the different 



estimates of a obtained by using the free and Nambu-Goto spectra might provide an 
estimate of systematic error on the determination of a from the lowest torelon mass 
due to our partial knowledge of the effective QCD string action. For l a = 3 one has 
a NG /a fs - 1 w 0.007. 



We expand the large-t behavior (|2l7| ) of G(t, 9) as 

G P (t, 9) « B e- Eot [1 + 9 2 h(t) + ...] (2.13) 

where we set 

B(9) = B + 9 2 B 2 + (2.14) 

E(9) = E + 9 2 E 2 + (2.15) 

and 

h(t) = ^-E 2 t (2.16) 

-DO 



Comparing the Eq. ( |2.13|) with the Eq. fl2.6|) , we find that 



Thus E 2 can be estimated from the difference 

Ah(t) = h{t)-h{t + a), (2.18) 

indeed 

lim Ah(t) = E 2 a (2.19) 

t— >oo 

Corrections are exponentially suppressed as exp[— (Eq — E )t] where Eq is the mass of 
the first excited state at 9 = 0. Assuming the free-string spectrum (|2.9|) , Eq — Eq = 
47i /L. Notice that, although Eq — E Q — > for L — > oo, this difference is not small in 
our calculations. Indeed, since we choose the lattice size L so that l a = y/aL ~ 3, 
(Eq - Eq)/Eq w Att/II w 1.4. 

Finally, the coefficient S2 of the 0(9 2 ) term in the expansion ( |1.6|) is obtained by 



S2 = ^r (2.20) 

s 2 is a dimensionless scaling quantity. It is expected to approach a constant in the 
continuum limit, with 0(a 2 ) scaling corrections. 

An analogous procedure can be used to determine the leading 0{9 2 ) term in the 
9 expansion of the lowest glueball mass M(9) around 9 = 0, cf. Eq. ( |1.7| ). In this 
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N 


(3 


lattice 


stat 


2 

a a 


a M Q ++ 




3 


5.9 


12 3 x 18 


25M/20 


0.0664(6) 


0.80(1) 


3.09(4) 


3 


6.0 


16 3 x 36 


25M/40 


0.0470(3) 


0.70(1) 


3.23(4) 


4 


10.85 


12 3 x 18 


16M/50 


0.0646(6) 


0.76(1) 


2.99(5) 


6 


24.5 


8 3 x 12 


9M/50 


0.114(2) 


0.83(1) 


2.46(4) 



Table 1: Some information on our Monte Carlo simulations for iV = 3, 4, 6. The estimates 
of the string tension a are obtained using Eq. ( |2.8| ). 



case we employ wall-wall correlators of plaquette-like operators with up to 6 links, all 
in spatial directions, in order to determine the ++ glueball mass. Correspondingly 
we define 

Ak(t) = k(t) - k(t + a) (2.21) 

where k(t) is the function analogous to h(t) defined from the glueball wall-wall cor- 
relators. Then, the 0(6 2 ) coefficient g% in the expansion ( |1.7|) can be obtained by 



92 = \, lim Ak(t) (2.22) 

where M is the ++ glueball mass. 

Finally, we mention that in order to improve the efficiency of the measurements 
we used smearing and blocking procedures (see, e.g., Refs. |5R El|) to construct new 
operators with a better overlap with the lightest propagating state. Our implemen- 



tation of smearing and blocking was already described in Ref. |26[ . We only mention 
that we constructed new super-links using three smearing, and a few (2-4) blocking 
steps, according to the value of L. These super-links were used to compute improved 
Polyakov lines or plaquette operators. 



3. Results 

In this section we present the results of our exploratory study of the 9 dependence 
of the spectrum using the method outlined in the preceding section. Table [I] con- 
tains some information on our MC runs for N = 3,4, 6 on lattices L 3 x T. Since 
the coefficients of the 9 expansions (|1.6[) and (|1.7|) are computed from connected 



correlation functions, such as ( |2.6|) , and turn out to be quite small, high statistics 
is required to distinguish their estimates from zero: Our runs range from 9 to 25 
million sweeps, with measures taken every 20-50 sweeps. This requirement repre- 
sents a serious limitation to the possibility of performing runs for large lattices and 
in the continuum limit, especially for large values of N, due also to the severe crit- 
ical slowing down discussed in Sec. |2.1| . For all values of (3 considered in this work, 
the autocorrelation time satisfies tq < 100 sweeps @. Furthermore, (3 values were 
chosen to lie in the weak-coupling region, i.e., beyond the first order phase transition 
in the case N = 6, and beyond the crossover region characterized by a peak of the 
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Figure 1: Plot of Ah(t) and Ak(t) for AT = 3 at (3 = 5.9 (above) and (3 = 6.0 (below). 
The data for Ak(t) are slightly shifted along the t axis. 



specific heat for iV = 3,4; see Ref. |26| for a more detailed discussion of this point. 
Runs generally started from cold configurations, to avoid problems due to metastable 
states at the transition (in the case N — 6). The lattice size L was chosen so that 
l a = tJoL > 3, which should be sufficiently large to obtain infinite-volume results 
(see, e.g., Refs. [|4], p6[]), at least within our precision. Due to the above-mentioned 
limitations, and in particular for N = 4, 6 we could afford only one value of (3, so 
that no stringent checks of scaling could be performed. For this reason our study 
should be still considered as a first exploratory investigation. 

Figs. [I] and ^| show the results for the discrete differences Ah(t) and Ak(t), cf. 
Eqs. (fTIED and ( ggg ), for N = 3 at (3 = 5.9,6.0 and for iV = 4, 6 respectively. As 
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expected, the signal degrades rapidly with increasing t. Anyway, they appear rather 
stable already for small values of t. As already discussed in Sec. |272|, the approach to 
a constant in the large-t limit should be exponential, as exp[— (Eq — E )t], where Eq 
is the energy of the first excited state at 9 = 0. In the case N = 3 the data at (3 = 6 
appear to approach the asymptotic behavior more rapidly than at (3 = 5.9. This 
should be due to the fact that a more effective blocking procedure can be applied 
when L = 16, rather than L = 12, achieving a better overlap with the lowest state. 

We estimate the coefficients s 2 and g 2 of the 0(8 2 ) terms in the expansions ( |1.6| ) 
and (|1.7| ) from the corresponding discrete differences Ah(t) and Ak(t) [cf. Eqs. (|2.19| ), 
( |2.20|) , and (|2.22|) 1, taking the data at t/a = 2 in the N = 3,4 runs, and at t/a = 1 



for iV = 6. In Table ^| we report the results. The estimates of s 2 and g 2 are small 
in all cases, and decrease with increasing N. For iV = 3 the results at (3 = 5.9 and 
(3 = 6.0 are consistent, supporting the expected scaling behavior. As final estimate 
one may consider 

s 2 = -0.08(1), g 2 = -0.06(2) for N = 3 (3.1) 
One may also consider the 9 dependence of the scaling ratio 

" + + ...), (3.2, 



where c 2 = g 2 — s 2 /2. Using the numbers reported in Table 0, we see that the 0(6 2 ) 
terms tend to cancel in the ratio. Indeed, we find c 2 = —0.02(2), —0.01(3), —0.01(2) 
respectively for iV = 3, 4, 6. 

For N = 4, 6 there are other independent A;-strings associated with representa- 
tions of higher n-ality. Analogously to the fundamental string, one may write 

a k (6) = a k (l + s kt2 6 2 + ...), (3.3) 

where a k is the fc-string tension at 9 = 0. The case k = 1 corresponds to the 
fundamental string tension, i.e. o\ = o and si >2 = s 2 . One may also consider the 
ratio R k = crk/cr, 

R k (6) = R k (l + r k , 2 8 2 + ...), (3.4) 



where R k is the ratio at 6 = (see e.g. Refs. [^3|, [26|, ^7], for recent numerical 
studies of the /c-string spectrum), and r kt2 = s k , 2 — s 2 . In the case iV = 4 there 
is one independent k string, <r 2 , besides the fundamental one; in the case N = 6 
there are two. Our results for the k > 1 strings are less stable. We obtained 
sufficiently precise results only from the simulation for N = 4. In the channel 
of Polyakov lines corresponding to the k = 2 string, we found a 2 a 2 = 0.091(2), 
Ah 2 (t/a = 1) = —0.044(13), while at distance t/a = 2 the signal was already 
unreliable. This leads to the estimate s 2 ,2 = —0.040(12). Note that s 2 ,2 ~ si,2 
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Figure 2: Plot of Ah(t) and Ak{t) for N = 4 at (5 = 10.85 (above) and N = 6 at = 24.5 
(below). The data for Ak(t) are slightly shifted along the t axis. 



N 





S2 


92 


3 


5.9 


-0.077(8) 


-0.05(2) 


3 


6.0 


-0.077(15) 


-0.07(4) 


4 


10.85 


-0.057(10) 


-0.04(3) 


6 


24.5 


-0.025(5) 


0.006(15) 



Table 2: Results for the coefficients S2 and g2, as derived from the discrete differences at 
distance t/a = 2 for N = 3, 4, and at t/a = 1 for iV = 6. 



{si t 2 = s 2 is reported in Table |2|), suggesting an even smaller 0(6 2 ) term in the ratio 
# 2 , i.e. |r 2 , 2 | < 0.02. 



- 10 - 



In conclusion, the above results show that the 0(9 2 ) terms in the expansion 
around 9 = of the spectrum of SU(iV) gauge theories are very small, especially 
when dimensionless ratios are considered. Moreover, they appear to decrease with 
increasing N, and the coefficients do not show evidence of convergence to a nonzero 
value. This is suggestive of a scenario in which the 9 dependence of the spectrum 
disappears in the large- N limit, at least for sufficiently small values of 9 around 
9 = 0. General large-iV scaling arguments applied to the Lagrangian ( |1 . 1| ) indicate 
9 = 9/N as the relevant Lagrangian parameter in the large- N limit [|J. In the case of 
the spectrum, this would imply that 0{9 2 ) coefficients should decrease as 1/N 2 . This 
is roughly verified by our results, taking also into account that they may be subject 
to scaling corrections, especially those at iV = 4, 6. For example, in the case of the 
string tension, s 2 ~ s 2i 2/N 2 with s 2 ,2 ~ —0.9. Of course, further investigations are 
required to put this scenario on a firmer ground. 

The hypothesis of a simple 9 dependence in the large- N limit may be extended 
to finite temperature, up to the first-order transition point. Recent studies [^, [30] 
have shown that in the large- N limit the topological properties remain substantially 
unchanged up to the first-order transition point. 



4. 9 dependence in the two-dimensional CP model 

Issues concerning the 9 dependence can also be discussed in two-dimensional CP 



N-l 



models [[31], [32], which are an interesting theoretical laboratory. Indeed they present 
several features that hold in QCD: asymptotic freedom, gauge invariance, existence of 
a confining potential between non gauge invariant states (that is eventually screened 
by the dynamical constituents), and non-trivial topological structure (instantons, 9 
vacua). Moreover, unlike SU(iV) gauge theories, a systematic 1/N expansion can be 



performed around the large- N saddle-point solution pi| , [32], p3 |. 

Analogously to four-dimensional SU(iV) gauge theories, one may add a 9 term 
to the Lagrangian, writing 

N 1 

Co = y D fl zD^z + i9— e^d^A,,, (4.1) 

where z is a ^-component complex scalar field subject to the constraint zz = 1, 
= izd^z is a composite gauge field, and = + iA^ is a covariant derivative. 
The topological charge density is q(x) = ^- e^ u d^A u . Then one may study the 9 
dependence of the ground state and other observables. In the following we discuss this 
issue within the 1/N expansion, performed keeping g fixed. Simple large- iV scaling 
arguments applied to the Lagrangian ( |4.1| ) indicate that the relevant 9 parameter in 
the large- iV limit should be 9 = 9/N. 
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As mass scale we consider the zero-momentum mass 1 M defined from the small- 
momentum behavior of the Fourier trasformed two-point correlation function of the 
operator Pij(x) = Zi(x)Zj(x), 

G P (x-y) = (TrP(x)P(y)), (4.2) 

i.e. from the relation 

G P (p)- 1 = Z-\M 2 + p 2 + 0(p 4 )] (4.3) 

where Z is a renormalization constant. 

Analogously to SU(iV) gauge theories, the ground state energy F{6) depends on 
9. One may define a scaling ground state energy f{9) and expand it around 8 = 0, 

f(0) = M- 2 [F(6) - F(0)] = ]fe 2 + b ^ e2n ) ( 4 ' 4 ) 



71=1 



where F{9) is defined as in Eq. ( |1.2p , M is the mass scale at 9 = 0, and C,bj are 
constants. C is the scaling ratio x/M 2 at 6 = 0, where x is the topological suscepti- 
bility, i.e. the two-point correlation function of the topological charge density at zero 
momentum. The correlation function of the topological charge density, and in par- 
ticular the topological susceptibility, has been computed within the 1/N expansion 
§§ |§. We have 

C = X /M 2 = ^ + 0(1/N 2 ) (4.5) 

The coefficients &2n are obtained from appropriate 2n-point correlation functions of 
the topological charge density operators at 9 = 0. For example 

b2 = (4.6) 
1 r (QVo-3((QVo) 2 ], (4.7) 

where Q = J d 2 xq{x) is the topological charge. 

We refer to Ref. for a discussion of the 1/N expansion within CP^ -1 models, 
its set up, and the list of the corresponding Feynman rules. In Fig. |3| we show the 
1/iV-expansion Feynman diagrams contributing to \i a ^ leading order. The analysis 
of the Feynman diagrams of the connected correlations necessary to compute bm 
shows that they are suppressed in the large- N limit, as 

b 2n = 0(1/N 2n ). (4.8) 



* 4= v 



lr This quantity is more suitable for a l/A^-expansion than the mass scale determined from the 
large-distance exponential decay of Gp(x), due to its analytical properties in 1/N [p3[. 
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Figure 3: Diagrams contributing to the four-point connected function %Ai cf. Eq. ([4.7|), 
within the 1/N expansion (the corresponding Feynman rules can be found in Ref. |33| ), 



This implies that the ground-state energy can be rewritten as 

f(9) = Nf(9 = 9/N), (4.9) 



71=1 



where C = NC and b 2n = N &2n are 0(1) in the large- N limit. Note the analogy 
with the expected 9 dependence of the ground-state energy in SU(iV) gauge theories, 
cf. Eq. (|1.5| ). The calculation of the coefficients of the leading large- N terms is rather 
cumbersome. Here, we only report the results obtained for the leading terms of b 2 
and 64 

t 27 1 1830 /Jin , 

b 2 = --=-, h = =-. 4.10 

5 7 

Within the 1/N expansion one may also study the dependence of the mass M 
on the parameter 9. We write 

M(9) = M (l + m 2 9 2 + ...) (4.11) 

The coefficient m>2 can be extracted from the connected correlations 

(Tr P(x)P(0)Q 2 ) e=0 - (Tr P(x)P(0)) e=0 (Q 2 ) e=Q (4.12) 

The analysis of its diagrams giving the corresponding 1/N expansion indicates that 
m 2 is suppressed as 

m 2 = 0(1/N 2 ) (4.13) 

This confirms the general large- N scaling arguments indicating 9 = 9/N as the 
relevant parameter in the large- N limit, as in the scenario put forward for the four- 
dimensional SU(iV) gauge theories. 
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